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Abstract 

Among various potential consequences of rarity is genetic erosion. Neutral 
genetic theory predicts that rare species will have lower genetic diversity than 
common species. To examine the association between genetic diversity and 
rarity, variation at eight DNA microsatellite markers was documented for 14 
Acropora species that display different patterns of distribution and abundance in 
the Indo-Pacific Ocean. Our results show that the relationship between rarity 
and genetic diversity is not a positive linear association because, contrary to 
expectations, some rare species are genetically diverse and some populations of 
common species are genetically depleted. Our data suggest that inbreeding is 
the most likely mechanism of genetic depletion in both rare and common cor- 
als, and that hybridization is the most likely explanation for higher than 
expected levels of genetic diversity in rare species. A significant hypothesis gen- 
erated from our study with direct conservation implications is that as a group, 
Acropora corals have lower genetic diversity at neutral microsatellite loci than 
may be expected from their taxonomic diversity, and this may suggest a height- 
ened susceptibility to environmental change. This hypothesis requires validation 
based on genetic diversity estimates derived from a large portion of the genome. 
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Introduction 

As a consequence of their small population size, neutral 
population genetics theory predicts that rare species will be 
genetically less diverse than common ones (Kimura 1983). 
In general, a positive linear relationship is expected 
between genetic diversity and population size (Wright 
1931), whereby as a species expands its population size, 
there are commensurate increases in genetic diversity. 
More specifically, for a neutral locus, the expected poly- 
morphism at mutation-drift equilibrium is proportional to 
the effective population size {N^ - the number of breeding 
individuals). Thus, in populations with large Ne, high levels 
of genetic variation are maintained, and this maximizes 
adaptive potential (Frankham et al. 2010). Furthermore, 
the variation in selective pressure between habitats within 



reefs leads to slightly different local adaptations within a 
population, and this facilitates higher productivity or sta- 
bility in the face of disturbance (Palumbi et al. 2008). 

The process whereby genetic diversity is lost in small 
populations is called genetic erosion (Vrijenhoek 1985), 
and this has been documented in populations of both 
plants and animals (Nevo et al. 1984; Elstrand and Elam 
1993; Baskauf et al. 1994; Frankham 1996). The causal 
factors of genetic erosion are mostly a combination of 
strong genetic drift through founder effects or bottle- 
necks, directional selection, clonality, and/or high levels 
of inbreeding (Kimura and Ohta 1971; Avise 1994; Willi 
et al. 2006; Frankham et al. 2010). Genetic erosion is 
problematic because it tends to reduce the fitness of 
individuals in a population. Hence, disturbance events, 
outbreaks of pathogens (Coltman et al. 1999), and other 
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stochastic events can force genetically depleted 
populations to extinction (Goodman 1987; Elstrand and 
Elam 1993; Fagen et al. 2002; Frankham et al. 2010). 

Under low population size, there is also an elevated risk 
that favorable alleles may be lost or that deleterious alleles 
will be fixed and both of these processes diminish the 
ability of individuals in a population to adapt to, or 
survive in, changing environments (Lande and Barrowc- 
lough 1987). Considering that genetic diversity can have 



important ecological consequences at the population, 
community, and ecosystem levels (Hughes et al. 2008), 
and especially for threatened species (Spielman et al. 
2004), it is important that population viability of rare 
species is examined (Palumbi 2003), and information 
about genetic diversity is made available for conservation 
decision making (van Oppen and Gates 2006). 

Among the 845 species of zooxanthellate scleractinian 
coral, published estimates of genetic diversity exist for only 



Table 1. Summary of population genetic data available for zooxanthellate scleractinian corals. 



Family 



Species 



Reference 



Pocilloporidae Seriatopora hystrix Ayre and Dufty (1994); Ayre and Hughes (2000, 2004); Maier et al. (2005); Underwood et al. 

(2007); van Oppen et al. (2008); Noreen et al. (2009); Bongaerts et al. (2010); Starger 
et al. (2010); van Oppen et al. (2011b). 
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Acroporidae 


Aaopora palmata 


Baums et al. (200Ba, 2006); Reyes and Schizas (2010); Palumbi et al. (2012) 


Acroporidae 


Acropora tenuis 


Marquez et al. (2002); Underwood et al.(2007); Nishikawa (2008); Underwood (2009) 


Acroporidae 


Aaopora valida 


Ayre and Hughes (2000, 2004) 


Faviidae 


Plesiastrea versipora 


Rodriguez-Lanetty and Hoegh-Guldberg (2002) 


Faviidae 


Favia fragum 


Goodbody-Gringley et al. 2010; 


Faviidae 


Goniastrea aspera 


Nishikawa and Sakai (2003); Nishikawa (2008) 


Faviidae 


Goniastrea australiensis 


Miller and Ayre (2008b) 


Faviidae 


Goniastrea favulus 


IVliller and Ayre (2008a) 


Faviidae 


Favia fragum 


Goodbody-Gringley et al. (2010) 


Faviidae 


Platygyra daedalea 


Miller and Ayre (2008a) 


Faviidae 


Platygyra sinensis 


Ng and Morton (2003) 


Faviidae 


Montastrea annularis 


Foster et al. (2007, 2012) 


Faviidae 


Montastrea cavernosa 


Goodbody-Gringley et al. (2012) 


Faviidae 


Montastrea faveolata 


Baums et al. (2010) 


Faviidae 


Diploria strigosa 


Atchison et al. (2008) 


Pectiniidae 


Mycedium elephantotus 


Yu et al. (1999); Dai et al. (2000) 


Fungiidae 


Fungia fungites 


Gilmour (2002) 


Fungiidae 


Heliofungia actiniiformis 


Knittweis et al. (2008) 


Dendropliylliidae 


Balanopliyllia europaea 


Goffredo et al. (2004) 


Siderastreidae 


Siderastrea stellata 


Neves et al. (2008) 


Siderastreidae 


Siderastrea radians 


Neves et al. (2008) 


Poritidae 


Porites lobata 


Polato et al. (2010) 


Astrocoeniidae 


Madracis decactis 


Atchison et al. (2008) 


Agariciidae 


Pavona gigantea 


Saavedra-Sotelo et al. (201 1) 
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4.6% of species (n = 39) (Table 1). These population 
genetic studies suggest that some high latitude populations 
of common coral species are vulnerable to genetic erosion 
(Ayre and Hughes 2004; Underwood et al. 2009), but others 
are not (Noreen et al. 2009). Until now, the level of genetic 
diversity in rare coral populations has only been examined 
among species restricted to the Atlantic Ocean (Baums et al. 
2005a, 2006, 2010; Foster et al. 2007, 2012; Atchison et al. 
2008; Neves et al. 2008; Reyes and Schizas 2010; Palumbi 
et al. 2012; Goodbody-Gringley et al. 2012); and one species 
(Pavona gigantea) restricted to the far Eastern Pacific Ocean 
(Saavedra-Sotelo et al. 2011). Thus, the genetic diversity 
and level of inbreeding in rare Indo-Pacific corals remain to 
be tested, these being the focus of this study. 

Acropora (staghorn corals) is the model group for this 
study because an extensive literature exists on the global 
ranges of Acropora species and we have abundance data 
that enable us to estimate means global census sizes. 
Acropora are extremely susceptible to coral bleaching, 
changes in water quality, disease, and predation (Marshall 
and Baird 2000; Bruno et al. 2007; Pearson 1981). Fur- 
thermore, because Acropora spp. are particularly impor- 
tant for reef formation, ecosystem function, and 



biodiversity, and 50% of species in the genus are listed in 
elevated categories of threat on the lUCN red list 
(Carpenter et al. 2008), the genetic implications of rarity 
in Acropora have direct conservation significance. 

This project is the first to perform a comparative analy- 
sis of genetic diversity over a large number of coral taxa. 
We examine the level of genetic diversity in 14 species of 
Acropora from the Indo-Pacific Ocean (nine rare and five 
common) encompassing 25 populations from 11 
geographic locations to obtain insights into their genetic 
diversity. We test the null hypothesis that rare species 
have lower genetic diversity than closely related common 
congeners, and we generate new a hypotheses pertaining 
to the susceptibility of Acropora corals to environmental 
change. 

Methods 

Samples of 14 species (nine rare, five common - Table 2) 
were collected from 11 locations across the Indo-Pacific 
(Fig. 1). Considering that "rarity" can apply not only to 
patterns of abundance but also to distribution (Brown 
1984; Gaston 1994), in this study, we examine the link 



Table 2. Summary of species, population sample sizes, and number of loci included in the final analysis. 

Number 

Species Population Geographic region Sample size of loci 



Acropora 


Orpheus Island 


Central GBR 


25 


7 


microphthalma 


Maldives 


North Indian Ocean 


12 


7 




Seychelles 


South Indian Ocean 


22 


7 




Kimbe Bay 


Papua New Guinea 


25 


7 


A. valida 


Orpheus Island 


Central GBR 


29 


7 




Heron Island 


Southern GBR 


26 


7 




Kimbe Bay 


Papua New Guinea 


20 


7 


A. austera 


Maldives 


Indian Ocean 


29 


6 




Arno Atoll 


North Central Pacific 


18 


5 




Majuro Atoll 


North Central Pacific 


24 


5 




Majuro - 20 branches 


Central Pacific 


20 


5 




from single colony 








A. millepora 


Ningaloo Reef 


Indian Ocean 


34 


8 




Orpheus Island 


Central GBR 


27 


8 


A. horrida 


Orpheus Island 


Central GBR 


27 


8 


A. papillare* 


Ningaloo Reef 


East Indian Ocean 


31 


7 




Orpheus Island 


Central GBR 


20 


8 




Okinawa - Japan 


North Pacific 


14 


8 


A. pichoni* 


Kimbe Bay 


Papua New Guinea 


6 


7 




Chuuk Lagoon 


Central West Pacific 


6 


7 


A. spathulata* 


Orpheus Island 


Central GBR 


28 


7 


A. kirstyae* 


Orpheus Island 


Central GBR 


27 


8 


A. tortuosa 


Rongelap Atoll 


North Central Pacific 


12 


7 


A. jacquelineae* 


Kimbe Bay 


Papua New Guinea 


20 


7 


A. kimbeensis* 


Kimbe Bay 


Papua New Guinea 


14 


7 


A. rongelapensis* 


Rongelap Atoll 


North Central Pacific 


12 


7 


A. walindii* 


Kimbe Bay 


Papua New Guinea 


14 


8 



Species marked with asterisk are rare. 
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Figure 1. Sampling locations for population genetic analysis. 



between genetic diversity and both estimated global 
census size and maximum global range size. 

To determine which species have a restricted global dis- 
tribution, the maximum global range of the 14 species 
included in this study was quantified using the World- 
Wide Acropora Database, which has 25,000 records based 
on over 30 years of collections (Wallace 1999). Longitudi- 
nal and latitudinal limits for each species were determined 
from the database, and the range was approximated as 
elliptical in shape with an area given by: Latitudinal 
Range/2 x Longitudinal Range/2 x Pi. Species were 
described as rare if their range is 1/lOth or less of the 
Acropora species with the largest global range {A. valida). 
Estimates of global census size were calculated according 
to Richards et al. (2008). For ease of interpretation, rare 
species are marked with an asterisk (*) throughout text 
(e.g., A. pichoni*). 

The population sample sizes of the 14 species included 
in this study range from 6 to 34 individuals (Table 2). It 
is important to note that conducting population genetic 
studies on rare species is challenging for a number of rea- 
sons, the principal one being that it is exceedingly diffi- 
cult to obtain sample sizes large enough to warrant 
interpretations to be made about population-level trends. 
For corals, the difficulty is further exacerbated by the 
remote nature of the locations where rare Acropora spe- 
cies occur, and the difficulty in identifying rare corals to 
the species level. Thus, for some of the rare species exam- 
ined in this study, local populations are so small that it is 
not feasible to obtain larger population samples. Hence, 
the small sample sizes and somewhat Hmited number of 



species examined prevent a rigorous test of the associa- 
tion between genetic diversity and rarity; however, this 
study provides a foundation from which the level of 
genetic diversity in rare corals can be further explored. 

All molecular samples examined in this project have 
matching skeletal voucher specimens that were identified 
by the author and verified by Dr. Garden Wallace. Small 
branches (2-5 cm) were collected from individual 
colonies and stored in absolute ethanol. To minimize 
sampling across multiple recruitment cohorts and asexu- 
ally derived clone mates, colony sizes and spacing were 
standardized (20-50 cm colony size, >20 m between colo- 
nies). DNA was extracted from approximately 20 mg of 
coral branch according to Underwood et al. (2009 - 
Appendix 1). Precipitated DNA was resuspended in 
100 /(L 0.1 mol/L Tris pH = 9 and stored at -20°C. 

Variation at nine variable tandem repeats (microsatel- 
lite markers) was documented using markers previously 
developed for Acropora (Baums et al. 2005b; van Oppen 
et al. 2007) (Table 3). Microsatellite polymerase chain 
reaction (PGR) products were initially examined using 
denaturing gel electrophoresis on the (Gorbett 
GelScan2000, Sydney, Australia). Microsatellite PGR prod- 
ucts were visualized using fluorescently labeled forward 
primers and unlabelled reverse primers. Once it was con- 
firmed via initial GelScan screening that the microsatel- 
lites would cross-amplify, genotyping was undertaken 
following the procedure described below. 

Micro satellites were pooled into three multiplex 
reactions (Table 4). Each PGR primer was labeled with a 
different fluorescent dye (TET, HEX, or FAM) and alleles 
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Table 3. Primer sequences. 



Locus name Primer sequence (5'-3') 



Ami!2_U02 


F - 


ACAAAATAACCCCTTCTACCT 




R - 


CTTC ATCTCTAC AGC C GATT 


Ami!2_006 


F - 


1 1 /™ A/^^TA AAA A A ^T^T^ A ^ A A 

CTTGACCTAAAAAACTGTCGTACAA 




R - 


/" 1 1 A 1 1 A^TA A A A A^^A^^A^A^A ATA A /" 1 1 1 

GTTATTACTAAAAAGGACGAGAGAATAAC 1 1 1 


Amil5_028 


F - 


/~ /^~r/^ A A A A A T — T/^ A A A A T/" 

G GTC G AAAA ATTG AAAAGTG 




R - 


AT/~ A A /^T^ /~T — 1 — 1 — T/^ A T/^ 

ATCACGAGTCCTTTTGACTG 


Amil2_022 


F- 


CTGTGGCCTTGTTAGATAGC 




R- 


AGATnGTGTTGTCCTGCTT 


Amil2_23 


F- 


GCAAGTGTTACTGCATCAAA 




R- 


TCATGATGCTTTACAGGTGA 


Amil2_007 


F - 


TAATGAGCAAACTCATTCATGG 




R- 


CTTTT CCAAGAGAAGTCAAGAA 


Amil2_010 


F- 


CAGCGAnAATATTTTAGAACAGTTTT 




R- 


CGTATAAACAAATTCCATGGTCTG 


Amil2_012 


F- 


TTTTAAAATGTGAAATGCATATGACA 




R- 


TCACCTGGGTCCCATTTCT 



Table 4. Multiplex reactions. 





Locus 


Repeal lype 


Label 


Multiplex 1 


Amil2_002 


(TG)10 


HEX 




Ami2_006 


(CA)4TA(CA)4 


FAM 




Amil5_028 


(TCACA)7TCAC 


TET 






(TCACA)4TCACTCACTCACA 




Multiplex 2 


Amil2_022 


(AC) 10 


TET 




Amil2_23 


(AG)7 


HEX 


Multiplex 3 


Amil2_007 


(TG)7AG 


TET 




Amil2_010 


TA(TG)1 1 


FAM 




Amil2_012 


GA(CA)6GA(CA)2 


HEX 



were scored as PGR product size in base pairs. Where 
more than two bands were observed in an individual, PGR 
products were cloned for subsequent sequencing to ensure 
that peaks were true alleles and did not represent nonspe- 
cific amplification. Gonditions for the PGR included using 
150-200 ng of DNA template and 5 //L 2x Qiagen Multi- 
plex PGR kit master mix in a 10 reaction in the pres- 
ence of 1 i)L of each primer and 3.25 yuL of H2O. PGR 
profile consisted of the initial denaturation step of 15 min 
followed by 35 cycles of 94° for 30 sec, 50° for 90 sec, and 
72° for 60 sec. The mix was incubated at 60°G for 30 min. 
Three microliters of the PGR product was electrophoresed 
in a 2% TAE-agarose gel in Ix TAE buffer to assess the 
yield. Successful products were then cleaned using the 
Sephadex resin in the Whatman Unifilter 800 system. One 
microliter of the purified PGR product was transferred to 
a skirted 96-well plate and sent for genotyping at the JGU 
Advanced Analytical Gentre. Fragment analysis was con- 
ducted on the Amersham MegaBase. To minimize 
genotyping errors, aU automated scorings of alleles were 
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checked manually, and rerunning the clean PGR product 
cleared uncertainties. 

In cases where over two alleles were detected in 
genotyping, the quaUty of genotyping results was cross- 
checked using standard cloning and sequencing tech- 
niques. Unlabelled microsatellite PGR products were 
cloned using the ligation kit, pGEM T easy (Promega, 
Sydney, Australia) (5 /^tL ligation buffer, 1 /xL pGEM-T 
Easy Vector, 3 /.(L PGR product, 1 pL DNA ligase) and 
incubated for 1-4 h at room temperature or overnight at 
4°C. The bacterial cells were transformed with a ligated 
vector using 60 /^L of NM522 competent cells. Gultures 
were spun in a benchtop centrifuge for 5 min at 
4000 rpm. The supernatant was removed and DNA was 
isolated using the plasmid isolation protocol in the RBG 
Hyfield Plasmid Mini Kit. The concentration of DNA was 
determined using a spectrometer and a minimum of 1 ^ug 
of purified DNA was dried and sent to Macrogen Inc. 
(vww.macrogen.com) for sequencing using SP6 and 
M13F vector primers. 

Analysis 

Microsatellite alleles were scored as a simple function of 
PGR product size. Genotypes for all loci were manually 
scored from electrophoretic data. Gonformity to the 
expectations of Hardy-Weinberg equilibrium (HWE) 
were established using a chi-square test (Miller and 
Benzie 1997) and significance values were adjusted with 
Benjamini-Hochberg (BY) correction for multiple 
comparisons (Benjamini and Yekutieli 2001; Narum 
2006) in GenAlex (Peakall and Smouse 2005). Genepop 
on the web (Raymond and Rousset 1995) was used to test 
for Unkage between loci under the following Markov 
Ghain parameters: 1000 dememorization, 100 batches, 
and 10,000 iterations per batch. Descriptive statistics, 
including proportion of polymorphic loci (P), number of 
alleles per locus (A), and observed and expected heterozy- 
gosity, were calculated to illustrate the distribution of 
genetic diversity within and between populations (Lewis 
and Zaykin 2001) - Nei's measure was used to correct for 
uneven sample size in heterozygosity estimates. Allele 
richness was calculated in Fstat v 2.9.3 (Goudet 2001), 
and this program was also used to correct for the uneven 
sample sizes among the populations examined. Allelic 
diversity and standard genetic distance were computed 
according to Nei (1987), and significance was corrected 
for multiple pairwise comparisons (Benjamini and Hoch- 
bergl995). 

The extent of inbreeding was summarized by the 
inbreeding coefficient, _Fis, in Fstat on the Web. This 
inbreeding coefficient assesses the effects of nonrandom 
mating vnthin subpopulations, as a measure of reductions 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



1871 



Rarity and Genetic Diversity in Staghorn Corals 



Z.T. Richards & M.J.H. van Oppen 



in the heterozygosity of individuals. The presence of null 
alleles (inconsistent ampUfication of alleles due to 
mutations in the primer binding region) was assessed in 
Microchecker v 2.2.3 (van Oosterhout et al. 2004). The 
probability of identity via sexual reproduction v^as then 
examined by calculating the proportion of unique 
multilocus genotypes (MLGs) at each site (Ng:N) (as per 
Underwood 2009). In situations where multilocus 
matches were identified within species, one individual 
from each pair was removed from subsequent analyses so 
that each unique genotype was represented only once. 

Statistical differences in genetic diversity and level of 
inbreeding among rare and common species were 
determined using the Kruskal-Wallace test implemented 
in SPSS 17 with a single outlier excluded (for further 
discussion see heterozygosity results for A. rongelapensis). 
The relationships between range size/census size and 
allelic richness/expected heterozygosity were initially 
examined with Pearson's Correlation Coefficient in R, 
and regression analysis was used to compare the goodness 
of fit (r^). Assumptions of Unearity, normality, and 
homogeneity of variances were assessed through 
examination of residuals and variables. The significance of 
linear {y = a + bx) and polynomial relationships 
iy = a + hx + cx^2) was also examined in multiple 
regression. 

Results 

A total of 531 individuals in 14 species of Acropora were 
genotyped (see Table 2). Thirty-eight percent of initial 
genotype runs either failed or had multiple peaks, so 
these samples were genotyped a second time to resolve 
peaks and 4% were genotyped three times. Overall, 10 
species had 100% polymorphic loci (Table 5); however, 
locus Apam3_166 did not amplify or amplified poorly in 
all samples, so it was removed from the analysis. 
Amil2_007 and Amil2_012 also provided mixed results 
and did not amplify in some populations of some species. 
For example, Amil2_007 did not amplify A. papillare* 
from Ningaloo Reef, but did amplify in A. papillare* from 
Orpheus Island and lapan; Amil2_012 did not amplify in 
any A. austera populations, but worked for all other 
species examined. Species showing <100% polymorphic 
loci were: A. papillare*, A. walindii*, A. valida, and 
A. austera. In A. valida, two populations contained 100% 
polymorphic loci, but the third (Heron Island) did not 
(57% of the loci were polymorphic). 



Clonality 

The population with the highest degree of clonality was 
the Majuro A. austera population where a single MLG 



Table 5. Total number of alleles screened across all loci (W), number 
and percentage of private alleles, number of locus pairs in linkage 
disequilibrium (LD), percentage polymorphic loci and MLG - identical 
multilocus genotypes, and probability of identity through sexual 
reproduction (Ng/N). 

% 

Polymorphic 



0.69 
1 

0.83 

1 

1 

0.93 
1 

1 
1 

0.94 
1 
1 
1 

0.84 

1 

1 

0.89 

0.96 

1 
1 



was repeated four times (Table 5) and the probability of 
asexual reproduction was 17%. The population with the 
highest occurrence of asexual reproduction was the Heron 
Island A. valida population where Ng/N was 69%. An 
additional five species displayed evidence of asexual 
reproduction including the two Ningaloo reef populations 
(A. papillare* 7% and A. millepora 6%), A. jacquelineae* 
from Kimbe Bay (16%); and A. kirstyae* and A. spathula- 
ta* from Orpheus Island (11% and 4%, respectively). 
The remaining populations of five species were fully 
sexually produced {A. horrida, A. tortuosa*, A. pichoni*, 
A. kimbeensis*, and A. rongelapensis*). 



Species 


Population 


N 


LD 


loci 


MLG 


Acropora 


Kimbe Bay 


51 


2 


100 


0 


micropthalma 


Seychelles 


53 


2 


100 


0 




Maldives 


35 


2 


100 


0 




Orpheus 


36 


2 


100 


1 




Island 










A. valida 


Orpheus 
Island 


41 


0 


100 


0 




Heron 


17 


0 


57.14 


3 




Island 












Kimbe Bay 


42 


0 


100 


0 


A. austera 


Majuro 


12 


0 


50 


4 




Arno 


14 


0 


66.7 


0 




Maldives 


29 


0 


100 


0 


A. papillare* 


Ningaloo 


41 


0 


87.5 


2 




Orpheus 


51 


0 


100 


0 




Island 












Japan 


32 


0 


87.5 


0 


A. millepora 


Orpheus 
Island 


48 


0 


100 


0 




Ningaloo 


62 


0 


100 


2 


A. pichoni* 


Kimbe Bay 


32 


0 


100 


0 




Truk 


34 


0 


100 


0 


A. horrida 


Orpheus 
Island 


43 


2 


100 


0 


A. jacquelineae* 


Kmibe Bay 


26 


0 


100 


3 


A. kimbeensis* 


Kimbe Bay 


37 


1 


100 


0 


A. tortuosa* 


Rongelap 
Atoll 


24 


2 


100 


0 


A. kirstyae* 


Orpheus 
Island 


49 


3 


100 


3 


A. spathulata* 


Orpheus 
Island 


28 


0 


100 


1 


A. walindii* 


Kimbe Bay 


18 


0 


62.4 


0 


A. rongelapensis* 


Rongelap 
Atoll 


36 


2 


100 


0 
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Figure 2. Chromatogram showing multiple peaks in locus Amil2_022 in Acropora valida from Orpheus Island. The examples of duplicated alleles 
that we describe could be evidence of hybridization events; however, duplication was restricted to a single locus (Amil2_022), and does not 
appear to represent a genome-wide (polyploidization) pattern. Duplication events cannot be explained as scoring errors or PCR artifacts because 
cloning and sequencing verified genotyping results. 



Potential polyploidy or multicopy loci 

Genotyping showed more than two peaks (three to five) 
in 15% of the individuals sampled (73 of the 531) 
(Fig. 2), suggesting either some of the loci are not 
single-copy in some species or that some species or popu- 
lations are polyploid. Species showing such patterns 
included A. microphthalma, A. valida, A. austera, A. 
kirstyae*, A. kimbeensis*, and A. pichoni*. Cloning and 
sequencing verified that > 2 alleles were present for locus 
Amil2_022 in A. valida and A. kimbeensis*, but not for 
any of the other loci or species. This suggests that this 
locus has undergone duplication in these species, rather 
than these species being polyploid. The existence of 
multiallelic profiles is not accommodated in computer 
programs that treat codominant markers, so this 
prevented the inclusion of locus Amil2_022 data in the 
analysis of the population within which they occur; 
however, it was included for the other species that 
showed < 2 alleles per individual at this locus. 

Heterozygosity 

Twenty-nine percent of samples displayed significantly 
lower observed heterozygosity than expected under HWE 
at P < 0.05. Two common species (A. millepora and 
A. valida) have the greatest proportion of loci with 
significant heterozygote deficits (number of loci in 
deficit = 62% and 70%, respectively). Significant hetero- 
zygote deficits were also detected in rare species (e.g., 
A. papillare*, A. pichoni*, A. kimbeensis*, and 
A. spathulata*) and null alleles were encountered 45 times 
(Appendix 1). 

Heterozygote deficits due to null alleles were corrected 
in 73% of cases. In one case, heterozygote deficit was due 
to large allele dropout, whereby shorter alleles are prefer- 
entially amplified, resulting in the less efficient amplifica- 
tion of large alleles; however, these data were corrected. 
For the remaining 27% of cases, significant deficits 
remained after correction for null alleles, suggesting that 
there are additional reasons for the deficits or there were 



not enough data to correct the null alleles. Correction 
reduced Fis scores and increased the number of popula- 
tions with heterozygote excess. For example, all three 
populations of A. papillare* changed to heterozygote 
excess after correction for null alleles (note: repeated 
MLG's were removed before analysis), while the Majuro 
A. austera population remained in deficit after null alleles 
were corrected. No null alleles were detected in A. ronge- 
lapensis* nor A. jacquelineae* . Significant heterozygote 
excess was detected in A. papillare* at 3 loci and in A. 
rongelapensis* at 6/7 loci with 100% observed heterozy- 
gosity recorded at 3 loci. 

Significant genotypic linkage disequilibrium was found 
in 18 of the performed tests (P < 0.05). Among common 
species, there was no strong link between loci with 
statistically significant genotypic linkage disequilibrium 
and populations in HWE disequilibrium. For example, 
two locus pairs (Amil2_002 and Amil5_002; Amil5_002 
and Amil2_010) were linked in all A. microphthalma 
populations; however, only one third of the 
A. microphthalma populations showed significant devia- 
tions from HWE (Appendix 1). However, for rare species, 
loci with statistically significant genotypic linkage disequi- 
librium also showed significant heterozygote deficits (e.g., 
A. tortuosa* and A. kirstyae*). The presence of linkage 
disequilibrium in association with heterozygote deficits 
may indicate inbreeding or it may be a sign that members 
of different populations have been sampled (i.e., Wahlund 
effect). 

Patterns of genetic diversity in rare and 
common species 

When both metrics of genetic diversity are plotted 
together with species ranked from most common to rare, 
it is obvious that the highest levels of expected heterozy- 
gosity occurred in rare species (Fig. 3). This figure also 
illustrates that among common species, the level of 
expected heterozygosity and allelic richness was similar; 
however, the results from these two genetic diversity met- 
rics were quite different for some of the rare species 
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Figure 3. Genetic diversity of rare and common Acropora species. Data were pooled for species sampled across multiple populations. Species are 
listed from most widespread (A. valida) to most geograplnically restricted {A. rongelapensis*). He, expected Ineterozygosity; AR, allelic riclnness. 



(Fig. 4). For example, in A. pichoni* mean allelic richness 
was particularly low, while expected heterozygosity was 
high. Generally patterns of allelic richness and expected 
heterozygosity at individual loci were extremely variable 
within species; hence, the large standard error bars. 

Acropora walindii* exhibited low allelic richness across 
all eight loci, suggesting that it is genetically eroded 
(Table 6). Furthermore, allelic fixation was detected at 
three loci for A. walindii*. Allele fixation was not 
restricted to rare species and additional examples of fixed 
loci are present in the Heron Island A. valida, Japan A. 
papillare*, and Majuro A. austera populations (Table 6). 
The most common locus to be fixed was Amil2_023. 
Mean allelic richness was greatest in A. microphthalma 
and this is driven largely by the 20 different alleles 
detected at locus Amil5_028 in the Seychelles population 
(Appendix 2). Mean expected heterozygosity was greatest 
in A. rongelapensis*, which was unexpected because it is 
the rarest of all species examined in this study. 

No significant difference (P < 0.05) could be detected 
in the level of expected heterozygosity {U = 16, 
P = 0.524) nor allelic richness (17=13.5, P = 0.230) 
between species that are numerically rare versus common. 
We tested whether there was a significant difference in 
the level of expected heterozygosity and allelic richness in 
species that are geographically restricted versus wide- 
spread, and again no significant difference (P < 0.05) was 
apparent (U = 19, P = 0.386 and U = 10.5, P = 0.109). 
We further examined the strength of linear and nonlinear 
relationships between genetic diversity and rarity metrics 
and the strength of the relationship (Table 7), and found 
that no significant linear or nonlinear relationships exist 



between expected heterozygosity and global range size or 
global census size (Figs. 4a and b). Similarly, there are no 
significant linear or nonlinear relationships between allelic 
richness and global range size or global census size 
(Figs. 4c and d). 

Inbreeding 

After correction for null alleles, Pjs values were extremely 
variable ranging from 0.54 to 0.215. The most highly 
inbred population was the Heron Island A. valida popula- 
tion followed by the Kimbe Bay A. pichoni* population 
(Fig. 5). There was no significant difference in the level of 
inbreeding between rare and common species {U = 21, 
P = 0.841) or between geographically restricted and wide- 
spread species (C7 = 19, P = 0.641). 

Discussion 

Our comparison of genetic diversity shows that, in 
general, there is a large amount of variability in the level 
of genetic diversity in Acropora spp with varying degrees 
of rarity. Contrary to expectations of neutral genetic the- 
ory, the rare Acropora species studied here do not have 
significantly lower genetic diversity (or higher levels of 
inbreeding) than common species. Thus, the association 
between the degree of rarity and genetic diversity in 
Indo-Pacific Acropora is not a simple positive linear 
relationship. For example, A. walindii* is restricted to 
Kimbe Bay in Papua New Guinea (Wallace 1999) and it 
has an estimated global population size of only 
1231 ± 615 individuals (Richards et al. 2008). A. walindii 
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Table 6. Patterns of allelic richness between species/populations across the eight loci examined. 



Species (min. pop. size) Population Amil2_02 Amil2_06 Amil5_028 Amil2_022 Amil2_023 Amil2_07 Amil2_010 Amil2_012 



Acropora pkhoni* (6) 


Kim be Bay 


2.879 


2.404 


3.564 


NA 


2.24 


2 


2.182 


3.309 




Chuuk 


3.309 


2.715 


3.343 


NA 


1.745 


3 


3.233 


3.309 




Overall mean 


3.12 


2.67 


3.499 


NA 


1.965 


2.414 


2.755 


3.208 


A. millepora (18) 


Palm Islands 


2.993 


3.908 


3.948 


13.305 


3 


3 


11.224 


1.91 




Ningaloo 


4.093 


4.296 


5.791 


9.118 


5.431 


5.937 


10.837 


6.292 




Overall mean 


3.575 


4.342 


6.074 


11.024 


4.431 


7.263 


13.845 


5.306 


A. austera (12) 


Majuro 


4 


4 


4 


NA 


1 


NA 


1 


NA 




Arno 


4.714 


3.857 


3.571 


NA 


1.999 


NA 


1 


NA 




Maldives 


5.659 


4.956 


2.807 


NA 


2.994 


NA 


5.039 


NA 




Overall mean 


5.772 


5.659 


4.434 


NA 


2.899 


NA 


4.311 


NA 


A. papillare* (14) 


Ningaloo 


2 


3.382 


5.453 


7.789 


6.228 


NA 


6.131 


4.444 




Palm Islands 


5.993 


2.7 


7.23 


11.13 


4.378 


NA 


6.231 


3.695 




Japan 


6 


4 


4 


4 


4 


NA 


5 


1 




Overall mean 


6.659 


6.307 


8.86 


10.781 


5.476 


NA 


7.162 


5.289 


A. microphthalma (19) 


Kim be Bay 


10.916 


4.706 


5.598 


NA 


5.892 


5.598 


8.921 


5.455 




Seychelles 


7.696 


7.68 


1 1 .045 


NA 


3 


6.711 


10.909 


5.846 




Maldives 


3 


6 


7 


NA 


4 


5 


6 


6 




Palm Islands 


5.711 


4.848 


3.711 


NA 


3.968 


4.848 


5.832 


5.711 




Overall mean 


10.8 


7.086 


9.288 


NA 


5.198 


7.554 


10.115 


6.708 


A. valida (20) 


Palm Islands 


4.287 


5.597 


12.168 


PP 


3.505 


NA 


2.69 


1.907 




Heron Island 


6.808 


3 


3 


PP 


1 


NA 


1.997 


1 




Kimbe Bay 


6 


6 


15 


PP 


4 


NA 


5 


2 




Overall mean 


7.022 


6.629 


13.843 


PP 


3.608 


NA 


4.621 


2 


A. kirstyae* (27) 


Palm Islands 


2 


3.994 


6.994 


13.65 


6.825 


5 


4.957 


5 


A. kimbeensis* (14) 


Kimbe Bay 


6 


4 


10 


PP 


3.995 


2 


6.786 


3.929 


A. horrida (27) 


Palm Islands 


2.617 


7.292 


8.347 


7.761 


3.478 


2 


3.943 


4.833 


A. walindii* (14) 


Kimbe Bay 


1.76 


1.791 


2.624 


2.002 


1 


1 


1.286 


1 


A. jacquelineae* (19) 


Kimbe Bay 


3 


2.636 


3.283 


NA 


4.895 


3 


4 


2 


A. spathulata* (26) 


Palm Islands 


1.462 


2.924 


3.993 


3.92 


2 


NA 


7.113 


4.33 


A. rongelapensis* (12) 


Rongelap Atoll 


6 


5 


8 


NA 


3 


4 


5 


5 


A. lorluosa* (12) 


Rongelap Atol 


4 


3 


4 


4 


2 


NA 


2 


5 



A wide range of variation in allelic richness was encountered within and between species. Allelic richness could not be calculated for some popu- 
lations at loci Amil2_022 because more than two alleles were encountered and we hypothesize that locus-specific duplication events have 
occurred. PP, polyploid; NA, not applicable. 



has exceedingly low genetic diversity and allele fixation 
is evident at 3 loci. Conversely, A. rongelapensis* has 
high genetic diversity with heterozygote excess at 6/7 loci 
and 100% observed heterozygosity (Hq) at three loci 
despite having an exceptionally small estimated global 
census size of only (224 ±117 individuals) (Richards 
et al. 2008). 

Our results supplement a growing number of rare/com- 
mon comparative studies in plants that challenge the 
traditional neutral genetic theory by demonstrating the 
absence of a significant positive relationship between rarity 
and genetic diversity (see Karron 1987, 1991; Hamrick and 
Godt 1990; Young and Brown 1996; Gitzendanner and Sol- 
tis 2000; Ellis et al. 2006). Factors which can lead to signifi- 
cant deviations from neutral theory in finite populations 
have recently been reviewed (Frankham 2012) and include 
balancing selection, selective sweeps and background selec- 
tion. While the practical limitations of this study precludes 



a robust test of neutral genetic theory, our results provide a 
novel attempt at integrating population demography and 
genetics, which is urgently needed to address issues related 
to the effects of climate change on species range and persis- 
tence (Lavergne et al. 2010). 

The wide range of variation we found between species 
with differing degrees of rarity undoubtedly reflects a 
range of population genetic factors. The fixed heterozy- 
gosity at three loci in A. rongelapensis*, the rarest species 
we examined in this study, may be explained in various 
ways. First, while asexual reproduction could result in 
fixed heterozygosity, no identical MLGs were identified. 
Second, high heterozygosity may reflect an old, stable, 
and persistent population or recruits that are derived 
from various genetically divergent sources (van Herwer- 
den et al. 2009). However, these explanations are not rele- 
vant here as A. rongelapensis* is a member of the large 
terminal dade in the Acropora phylogeny that is consid- 
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Table 7. Regression statistics showing the strength and significance 
of linear and non-linear relationships between range size or census 
size and allelic richness or expected heterozygosity. 













Sig. 


Association 


Relationship 




df 


P 


<0.05 


Range size and 


Linear 


0.211 


12 


0.098 


NS 


allelic 


Polynomial 


0.235 


11 


0.229 


NS 


richness 








(x = 0.383; 
x^ = 0.570) 




Range size 


Linear 


0.015 


12 


0.680 


NS 


and expected 


Polynomial 


0.235 


11 


0.229 


NS 










(x = 0 1 1 7' 
x^ = 0.103) 




Census size and 


Linear 


0.159 


12 


0.159 


NS 


allelic richness 


Polynomial 


0.159 


11 


0.386 
(x = 0.809; 
x^ = 0.978) 


NS 


Census size and 


Linear 


0.005 


12 


0.812 


NS 


expected 


Polynomial 


0.056 


11 


0.729 


NS 


heterozygosity 








(x = 0.480; 
x^ = 0.457) 





NS, nonsignificant. 



ered relatively young (<5 my - Richards et al. 2008) and 
only a small number of isolated populations of A. 
rongelapensis* have been located across its continuous 
distribution range. 

The preferred explanation for the finding of 100% 
observed heterozygosity at 3 loci in A. rongelapensis* is 
interspecific hybridization. Hybridization has been dem- 
onstrated previously in Acropora communities, namely in 
the Caribbean where only three extant species exist, and 
A. prolifera was found to be the product of hybridization 
between A. palmata and A. cervicornis (van Oppen et al. 
2000; VoUmer and Palumbi 2002). We hypothesize that 
the individuals of A. rongelapensis examined in this study 
are Fl hybrids; however, the parental lineages have not 
been established to date. It is likely that hybridization 
may also have contributed to the higher than expected 
genetic diversity in other rare species (e.g., A. kimbeen- 
sis*). The proposition that some rare species have hybrid 
ancestries is further supported by a phylogenetic analysis 
of nuclear and mitochondrial DNA that shows that some 
rare species are monophylotic for mitochondrial DNA 
and polyphyletic for nuclear DNA (Richards et al. 2008). 

Furthermore, the finding of gene duplication events at 
Amil2_022 in A. kimbeensis* supports the suggestion that 
hybridization is a mechanism driving the evolution of 
genetic diversity in rare species. Single locus duplication 
events have been reported in another microsateUite study 
involving A. millepora (Wang et al. 2008). Chimerism 
(where juveniles settle together and fuse - Barki et al. 
2002; Puill-Stephan et al. 2009) and the retention of a 
polar body during fertilization (Baums et al. 2005b) were 



suggested as possible explanations for the occurrence of 
more than two alleles per locus. Chimerism is not likely 
to explain our observations because DNA was examined 
from a small portion of a single branch only. It is possible 
that Amil2_022 is not a single copy marker; however, the 
finding that duplication was restricted to a few individu- 
als indicates that this may also reflect a recent region-spe- 
cific dupUcation event caused by transposable element 
activity, replication slippage, or aberrant crossing over 
(Bennetzen 2002). Another possible mechanism for 
tri-aUelic patterns is somatic mutations that result in 
genetic mosaics (van Oppen et al. 2011a). We suggest 
that gene/genome duplication events should be further 
examined as mechanisms that drive the evolution of 
genetic diversity in corals (for further discussion of gene 
duplication see Stebbins 1940; Ohno 1970; Lynch and 
Conery 2000; Zhang 2003). 

Populations with low genetic diversity (such as A. 
walindii*, A. pichoni*, and A. jacquelineae* in Kimbe Bay 
and A. papillare* in Japan) may be more vulnerable and 
have a higher probability of local extinction following dis- 
turbance. Local extinction (the disappearance of a species 
from part of it range) and ecological extinction (when a 
species is reduced to such low abundance that, although 
still present, it no longer plays its typical ecological role) 
are precursors to global extinction. Data here suggest that 
some populations of common Acropora species may be 
more vulnerable to local extinction events than previously 
thought. For example, the A. valida population at Heron 
Island had low genetic diversity, significant heterozygote 
deficits, and some alleles are fixed. Heterozygote deficits 
are commonly caused by clonality, nuU alleles, inadver- 
tently sampling disparate populations (i.e., Wahlund 
effects, Wahlund 1928) or inbreeding (Wright 1922). In 
this analysis, clonality and nuU alleles can be excluded as 
possible explanations for the observed deficits. Moreover, 
evidence presented suggests that inbreeding is the primary 
explanation for heterozygote deficits (e.g., for A. valida at 
Heron Island); however, we cannot exclude the possibility 
that Wahlund effects may have contributed to the 
observed deficits. 

The concept of corals having low genetic diversity and 
low A/e is counterintuitive, given that corals such as Acro- 
pora have the potential for very high fecundity, in addi- 
tion to high levels of outbreeding and gene flow (both of 
which renew genetic variation and increase N^, CabaUero 
and Hill 1992). A meta-analysis which tested for devia- 
tions from the predicted positive linear association 
between genetic diversity and population size suggests 
that species with high fecundity (such as fish, oysters, 
shrimp, and seaweed) have significantly reduced NJN 
ratios (Frankham 2012). Thus, despite most corals having 
high fecundity, they tend to experience very high mortali- 
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Figure 4. Correlation between range size/estimated global census size and mean expected heterozygosity/allelic riclnness for 14 species witin 
different distribution/abundance patterns. (A) A nonsignificant nonlinear relationslnip exists between maximum range size and mean He 
(r" = 0.235, df = 11, P-value X = 0.103). (B) A nonsignificant negative linear relationship exists between estimated global census size and mean 
He (r^ = 0.005; df = 12; P = 0.812). (C) A nonsignificant linear relationship exists between maximum range and allelic richness (f^ = 0.211, 
df = 12, P = 0.098). (D) A weak nonsignificant positive linear relationship exists between estimated global census size and mean allelic richness 
(r^ = 0.159, df = 12, P= 0.159). 



ties in the early life stages, and for any given year, most 
of the recruited young may be derived from a few large 
parents, hence there is high variance in progeny number, 
and consequently, is likely to be small (Hughes et al. 
1992; Hedgecock 1994). 

An unexpected outcome of this study was our finding 
that despite being a taxonomically diverse group in the 
Indo-Pacific, our overall mean estimate of genetic 
diversity is quite low in comparison with the mean non- 
Acropora scleractinian diversity estimate reported by 
Shearer et al. (2009) (7.94 alleles per species specific 
locus, see Fig. 6). After correction for uneven sample 



sizes, the mean number of alleles per locus reported 
across the 25 Acropora populations was 4.64 (±0.2 SE) 
alleles per locus. This level of genetic diversity is just over 
half of the "conservative mean" presented for non-Acro- 
pora corals and more closely resembles those reported 
from isolated high latitude populations (Miller and Ayre 
2008b). All of the mean estimates of genetic diversity pre- 
sented here for Indo-Pacific Acropora are well short of 
the number presented for A. palmata, a rare species 
endemic to the Atlantic Ocean (mean 14.4). 

To confirm the hypothesis that Acropora have lower 
genetic diversity than other types of corals, genetic 
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diversity estimates could be corrected for differential 
mutation rates, or genome scans based on, for example, 
single nucleotide polymorphisms could provide a more 
accurate picture of functional genetic diversity across a 
large portion of the genome. 

Field observations and mesocosm results in other 
systems jointly support the hypothesis that genotypic 
diversity is important in providing different responses to 
environmental variation (Whitham et al. 2003; Gamfeldt 
et al. 2005; Hughes and Stachowicz 2009; Wang et al. 
2012). However, while our results suggest that some pop- 
ulations of rare and common species have low genetic 
diversity, it would be erroneous to directly infer anything 
about the adaptive potential of these populations or the 
way they will respond to climate change. Microsatellite 
markers are generally considered neutral (Estoup and 
Angers 1998) and hence they are not functionally con- 
strained or under selection (Frankham et al. 2010). Thus, 
for the loci examined here, genetic drift and gene flow, 
rather than selection, determine their allele number and 
heterozygosity. Therefore, the finding of low genetic 
diversity at the relatively small number of neutral loci 
examined here may not be functionally informative. We 
recommend that further robust species- specific popula- 
tion estimates at neutral and functional loci are required 
to fuUy understand the relationship between genetic 
diversity, adaptive potential, and persistence. 

The conservation impHcations of this study are found 
in both theory and practice. Theoretically, we show that 
there is a large amount of variabiHty in the genetic diver- 
sity of Acropora corals, which appears to be driven by 
processes other than rarity. We do not go so far as to 
suggest that rarity has a minor impact on genetic diversity 
(see Gillespie 2001; Bazin et al. 2006); however, the find- 
ing of nondepletion in some rare species and depletion in 
some common species suggests that genetic diversity is 
governed by a complex range of factors. Practically, 
we emphasize that three of the species which are shovm 
here to have exceptionally low levels of genetic diversity 
(A. walindii*, A. jacquelineae* , and A. papillare* at Japan) 
are listed by the lUCN as "Vulnerable," meaning that 
they have an elevated risk of extinction risk this century 
(see Carpenter et al. 2008). The finding of low genetic 
diversity in these species further exacerbates their 
threatened status and we propose that targeted species 
monitoring and management intervention may be 
necessary to safeguard these species. 
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Appendix 1 . species, population, locus, and number of samples (N), alleles (A), expected (H^), and observed (Hq) 
heterozygosity. Inbreeding coefficient (fis), nuU alleles (NuU), corrected inbreeding coefficient (fisC). Large allele drop- 
out (LAD), not enough data (NED), not applicable (NA). Asterisks indicate significant deviations fi-om Hardy-Weinberg 

equilibrium (P < 0.05). 

Locus W A He Ho F,s Null F,sC 



Acropora microphthalma Kimbe Bay 2 23 

6 23 
28 23 

22 NA 

23 23 

7 22 
10 25 
12 25 

Seychelles 2 22 

6 22 
28 22 

22 NA 

23 22 

7 22 
10 22 
12 22 

Maldives 2 19 

6 19 
28 19 

22 NA 

23 19 

7 19 
10 19 
12 19 

Orpheus Island 2 22 

6 22 
28 22 

22 NA 

23 22 

7 22 
10 22 
12 22 

A. valida Orpheus Island 2 28 

6 25 
28 26 

22 NA 

23 27 

7 28 
10 28 
12 28 

Heron Island 2 17 

6 13 
28 14 

22 NA 

23 16 

7 NA 
10 17 
12 17 

Kimbe Bay 2 19 

6 19 

28 19 



12 


0.849 


0.826 


0.105 


No 


0.105 


5 


0.691 


0.783 


-0.147 


No 


-0.147 


6 


0.509 


0.435 


0.096 


No 


0.096 


NA 


NA 


NA 


NA 


No 


NA 


6 


0.778 


0.652 


0.136* 


No 


0.136 


6 


0.478 


0.045 


0.92* 


Yes 


0.119 


10 


0.728 


0.72 


0.031 


No 


0.031 


6 


0.406 


0.32 


0.232* 


No 


0.232 


8 


0.747 


0.611 


0.141 


No 


0.141 


7 


0.638 


0.529 


0.191 


No 


0.191 


20 


0.787 


0.389 


0.515* 


Yes 


0.169 


NA 


NA 


NA 


NA 


NA 


NA 


3 


0.575 


0.714 


-0.171 


No 


-0.171 


7 


0.702 


0.591 


0.181* 


No 


0.181 


12 


0.778 


0.545 


0.32* 


Yes 


0.267 


6 


0.635 


0.318 


0.516* 


Yes 


0.115 


3 


0.661 


0.579 


0.15 


No 


0.15 


6 


0.726 


0.789 


-0.061 


No 


-0.061 


7 


0.74 


0.684 


0.102 


No 


0.068 


NA 


NA 


NA 


NA 


NA 


NA 


4 


0.35 


0.421 


-0.176 


No 


-0.176 


5 


0.56 


0.333 


0.566 


Yes 


0.186 


5 


0.738 


0.533 


0.176* 


Yes 


0.16 


5 


0.616 


0.867 


-0.397 


No 


-0.397 


6 


0.624 


0.773 


-0.216 


No 


-0.216 


5 


0.636 


0.818 


-0.264 


No 


-0.264 


4 


0.17 


0.091 


0.485 


Yes 


-0.05 


NA 


NA 


NA 


NA 


NA 


NA 


4 


0.458 


0.273 


0.423* 


Yes 


0.133 


5 


0.522 


0.409 


0.238 


No 


0.238 


6 


0.667 


0.864 


-0.273 


No 


-0.273 


6 


0.674 


0.909 


-0.329* 


No 


-0.329 


6 


0.596 


0.357 


0.281* 


Yes 


0.028 


6 


0.728 


0.52 


0.238* 


Yes 


0.105 


13 


0.877 


0.654 


0.231* 


Yes 


0.225 


NA 


NA 


NA 


NA 


NA 


NA 


4 


0.141 


0.074 


0.378 


Yes 


0.043 


7 


0.694 


0.179 


0* 


No 


0 


3 


0.364 


0.321 


0.143 


No 


0.143 


2 


0.069 


0 


1* 


No 


1 


7 


0.78 


0.588 


0.223 


NED 


0.223 


3 


0.506 


0.077 


0.727* 


NED 


0.727 


3 


0.612 


0 


0.87* 


NED 


0.87 


NA 


NA 


NA 


NA 


NA 


NA 


1 


0 


0 


0 


NED 


0 


NA 


NA 


NA 


NA 


NED 


NA 


2 


0.111 


0.118 


-0.045 


NED 


-0.045 


1 


0 


0 


0 


NED 


0 


6 


0.735 


0.737 


0 


No 


0 


6 


0.719 


0.368 


0.535* 


Yes 


0 


15 


0.904 


0.421 


0.571* 


Yes 


0.085 
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A. austera 



Majuro 



Arno 



Maldives 



A. papillare 



Ningaloo 



Orpheus Island 



Japan 



A. millepora 



Orpheus Island 



Locus 


N 


A 


He 


Ho 




Null 


FisC 


22 


NA 


NA 


NA 


NA 


NA 


No 


NA 


23 


19 


4 


0.389 


0.368 


0.089 


No 


0.089 


7 


19 


4 


0.609 


0.105 


0* 


No 


0 


10 


19 


5 


0.503 


0.368 


0.3* 


No 


0.3 


12 


19 


2 


0.188 


0 




Yes 


0.3 


2 


12 


4 


0.462 


0.417 


0.141 


No 


0.141 


6 


12 


4 


0.66 


0.417 


0.405 


No 


0.405 


28 


12 


4 


0.601 


0.25 


0.612* 


Yes 


0.12 


22 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


23 


11 


1 


0 


0 


0 


No 


0 


7 


NA 


NA 


NA 


NA 


NA 


No 


NA 


10 


12 


1 


0 


0 


0 


No 


0 


12 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


2 


14 


5 


0.666 


0.857 


-0.253 


No 


-0.253 


6 


14 


4 


0.666 


0.857 


-0.253 


No 


-0.253 


28 


10 


4 


0.27 


0.3 


-0.04 


No 


-0.04 


22 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


23 


14 


2 


0.191 


0.214 


-0.083 


No 


-0.083 


7 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


10 


14 


1 


0 


0 


0 


No 


NA 


12 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


2 


28 


7 


0.733 


0.786 


-0.054 


No 


-0.054 


6 


28 


6 


0.625 


0.536 


0.161 


No 


0.161 


28 


28 


3 


0.307 


0.286 


0.087 


No 


0.087 


22 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


23 


29 


3 


0.511 


0.586 


-0.131 


No 


-0.131 


7 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


10 


29 


6 


0.728 


0.655 


0.118* 


No 


0.118 


12 


29 


6 


0.702 


0.448 


0* 


No 


NA 


2 


28 


3 


0.07 


0.036 


0.5* 


Yes 


-0.038 


6 


28 


4 


0.563 


0.179 


0.692* 


Yes 


-0.041 


28 


28 


6 


0.735 


0.357 


0.528* 


Yes 


0.028 


22 


25 


9 


0.788 


1 


-0.25* 


No 


-0.25 


23 


25 


7 


0.626 


0.6 


0.061 


No 


0.061 


7 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


10 


28 


7 


0.721 


0.571 


0.225* 


Yes 


0.015 


12 


26 


5 


0.57 


0.462 


0.208* 


No 


0.208 


2 


20 


7 


0.639 


0.75 


0.164 


No 


0.164 


6 


20 


3 


0.411 


0.2 


0.532* 


Yes 


-0.282 


28 


20 


8 


0.8 


0.8 


0.026 


No 


0.026 


22 


20 


14 


0.814 


1 


-0.204 


No 


-0.204 


23 


20 


5 


0.585 


0.15 


0.755* 


Yes 


-0.078 


7 


20 


3 


0.421 


0.15 


0* 


No 


0 


10 


20 


7 


0.706 


0.75 


-0.036 


No 


-0.036 


12 


20 


4 


0.554 


0.8 


-0.424 


No 


-0.424 


2 


14 


6 


0.755 


0.643 


0.185* 


No 


0.185 


6 


14 


4 


0.543 


0.571 


-0.015 


No 


-0.015 


28 


14 


4 


0.707 


0.071 


0.906* 


NED 


0.906 


22 


14 


4 


0.612 


0.857 


-0.368 


No 


-0.368 


23 


14 


4 


0.199 


0.143 


0.316 


No 


0.316 


7 


14 


4 


0.487 


0.643 


0 


No 


0 


10 


14 


5 


0.691 


0.714 


0.004 


No 


0.004 


12 


14 


1 


0 


0 


0 


No 


0 


2 


26 


3 


0.556 


0.385 


0.326 


No 


0.326 


6 


26 


4 


0.499 


0.385 


0.248* 


No 


0.248 
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Appendix 1. (Continued) 







Locus 


W 


A 


He 


Ho 


Fis 


Null 








28 


26 


4 


0.362 


0.192 


0.485* 


Yes 


0.162 






22 


26 


16 


0.873 


0.846 


0.05 


No 


0.05 






23 


26 


3 


0.473 


0.308 


0.367* 


Yes 


0.006 






7 


18 


3 


0.642 


0 


1* 


Yes 


1 






10 


26 


13 


0.749 


0.731 


0.044 


No 


0.044 






12 


26 


2 


0.074 


0.077 


-0.02 


No 


-0.02 




Ningaloo 


2 


32 


5 


0.409 


0.469 


-0.131 


No 


-0.131 






6 


32 


5 


0.445 


0.156 


0.658* 


Yes 


0.168 






28 


33 


6 


0.786 


0.394 


0.51* 


Yes 


0.155 






22 


33 


12 


0.805 


0.727 


0.112 


No 


0.112 






23 


33 


7 


0.597 


0.515 


0.152 


No 


0.152 






7 


32 


7 


0.754 


0.219 


0.718* 


Yes 


-0.084 






10 


33 


13 


0.865 


0.576 


0.348* 


Yes 


0.102 






12 


33 


7 


0.648 


0.333 


0.498* 


Yes 


0.022 


A. pichoni 


Kimbe Bay 


2 


6 


6 


0.694 


0.5 


0.362 


NED 


0.362 






6 


6 


3 


0.611 


0 


1* 


NED 


1 






28 


6 


8 


0.847 


0.833 


0.107 


NED 


0.107 






22 


NA 


NA 


NA 


NA 


NA 


NED 


NA 






23 


6 


4 


0.514 


0.5 


0.118 


NED 


0.118 






7 


2 


2 


0.375 


0.5 


0 


NED 


0 






10 


6 


3 


0.542 


0.167 


0.737 


NED 


0.737 






12 


6 


6 


0.806 


0.667 


0.259 


NED 


0.259 




Truk 


2 


6 


6 


0.806 


1 


-0.154 


NED 


-0.154 






6 


6 


5 


0.667 


0.5 


0.333 


NED 


0.333 






28 


6 


7 


0.806 


0.833 


0.057 


NED 


0.057 






22 


NA 


NA 


NA 


NA 


NA 


NED 


NA 






23 


6 


2 


0.375 


0.5 


-0.25 


NED 


-0.25 






7 


2 


3 


0.625 


0.5 


0.5 


NED 


0.5 






10 


5 


5 


0.78 


0.6 


0.333 


NED 


0.333 






12 


6 


6 


0.806 


0.333 


0.643* 


NED 


0.643 


A. horrida 


Orpheus Island 


2 


26 


3 


0.144 


0.038 


0.742* 


Yes 


0.059 






6 


26 


8 


0.803 


0.346 


0.582* 


Yes 


0.019 






28 


26 


10 


0.686 


0.346 


0.51* 


Yes 


0.067 






22 


22 


8 


0.834 


0.364 


0.579* 


Yes 


0.012 






23 


23 


4 


0.427 


0.565 


-0.303 


No 


-0.303 






7 


17 


2 


0.457 


0 


1* 


Yes 


-0.333 






10 


18 


4 


0.335 


0.389 


-0.133 


No 


-0.133 






12 


18 


5 


0.466 


0.611 


-0.285 


No 


-0.285 


A. jacquelineae 


Orpheus Island 


2 


17 


3 


0.631 


0.529 


0.191 


No 


0.191 






6 


17 


3 


0.258 


0.176 


0.342 


No 


0.342 






28 


17 


4 


0.306 


0.176 


0.448 


No 


0.448 






22 


NA 


NA 


NA 


NA 


NA 


No 


NA 






23 


19 


7 


0.578 


0.789 


-0.343 


No 


-0.343 






7 


11 


3 


0.657 


0.545 


0.216* 


No 


0.216 






10 


11 


4 


0.442 


0.455 


0.02 


No 


0.02 






12 


11 


2 


0.236 


0.091 


0.643 


No 


0.643 


A. kimbeensis 


Kimbe Bay 


2 


13 


6 


0.719 


0.385 


0.496* 


Yes 


0.163 






6 


13 


4 


0.666 


0.538 


0.229 


No 


0.229 






28 


13 


10 


0.858 


0.615 


0.319 


Yes 


0.319 






22 


NA 


NA 


NA 


NA 


NA 


No 


NA 






23 


14 


4 


0.459 


0.143 


0.708* 


Yes 


0.071 






7 


13 


2 


0.497 


0.462 


0.111 


No 


0.111 






10 


14 


7 


0.694 


0.857 


-0.2 


No 


-0.2 






12 


14 


4 


0.656 


0.429 


0.378 


No 


0.378 


A. tortuosa 


Rongelap Atoll 


2 


12 


4 


0.462 


0.25 


0.492* 


Yes 


0.135 
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Locus 


N 


A 


He 


Ho 




Null 








6 


12 


3 


0.559 


0.25 


0.582* 


Yes 


-0.021 






28 


12 


4 


0.608 


0.25 


0.616* 


LAD 


0.616 






22 


12 


4 


0.563 


0.917 


-0.603 


No 


-0.603 






23 


12 


2 


0.153 


0.167 


-0.048 


No 


-0.048 






7 


NA 


NA 


NA 


NA 


NA 


NA 


NA 






10 


12 


2 


0.278 


0.333 


-0.158 


No 


-0.158 






12 


12 


B 


0.778 


0.833 


-0.028 


No 


-0.028 


A. kirstyae 


Orpheus Island 


2 


24 


2 


0.499 


0.125 


0.759* 


Yes 


-0.004 






6 


24 


4 


0.582 


0.167 


0.724* 


Yes 


-0.129 






28 


24 


7 


0.777 


0.583 


0.269 


LAD 


0.269 






22 


23 


14 


0.765 


0.913 


-0.173 


No 


-0.173 






23 


23 


7 


0.433 


0.348 


0.218 


No 


0.218 






7 


22 


5 


0.249 


0.136 


0.471* 


Yes 


-0.122 






10 


23 


5 


0.631 


0.826 


-0.288 


No 


-0.288 






12 


23 


5 


0.632 


0.609 


0.06* 


No 


0.06 


A. spathulata 


Orpheus Island 


2 


26 


2 


0.038 


0.038 


0* 


No 


0 






6 


26 


3 


0.443 


0.077 


0.832* 


Yes 


-0.098 






28 


24 


4 


0.666 


0.5 


0.269* 


No 


0.269 






22 


13 


4 


0.388 


0.154 


0.628* 


Yes 


0.13 






23 


12 


2 


0.153 


0.167 


-0.048 


No 


-0.048 






7 


NA 


NA 


NA 


NA 


NA 


No 


NA 






10 


20 


8 


0.8 


0.9 


-0.1 


No 


-0.1 






12 


22 


5 


0.546 


0.636 


-0.142 


No 


-0.142 


A. walindii 


Kimbe Bay 


2 


14 


2 


0.408 


0.571 


-0.368 


NED 


-0.368 






6 


13 


2 


0.426 


0.154 


0.662 


NED 


0.662 






28 


12 


5 


0.681 


0.5 


0.305 


NED 


0.305 






22 


14 


3 


0.518 


0.5 


0.071 


NED 


0.071 






23 


14 


1 


0 


0 


0 


NED 


0 






7 


2 


1 


0 


0 


0 


NED 


0 






10 


14 


3 


0.135 


0.143 


-0.02 


NED 


-0.02 






12 


14 


1 


0 


0 


0 


NED 


0 


A. rongelapensis 


Rongelap Atoll 


2 


12 


6 


0.722 


1 


-0.347* 


No 


-0.347 






6 


12 


5 


0.726 


0.75 


0.01 


No 


0.01 






28 


12 


8 


0.83 


0.833 


0.039* 


No 


0.039 






22 


NA 


NA 


NA 


NA 


NA 


No 


NA 






23 


12 


3 


0.538 


1 


-0.846* 


No 


-0.846 






7 


12 


4 


0.726 


0.917 


-0.222* 


No 


-0.222 






10 


12 


5 


0.694 


0.667 


0.083 


No 


0.083 






12 


12 


5 


0.674 


1 


-0.451 


No 


-0.451 
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Appendix 2. Number of alleles per locus. 



Species 


Amil2_02 


Amil2_06 


Amil5_028 


Amil2_022 


Amil2_023 


Amil2_07 


Amil2_010 


Amil2. 


Acropora micropthalma KB 


12 


5 


6 


0 


6 


6 


10 


6 


A. micropthalma SY 


8 


7 


20 


0 


3 


7 


12 


6 


A. micropthalma ML 


3 


6 


7 


0 


4 


5 


5 


5 


A. micropthalma PI 


6 


5 


4 


0 


4 


5 


6 


6 


A. valida PI 


6 


6 


13 


0 


4 


7 


3 


2 


A. valida HI 


7 


3 


3 


0 


1 


0 


2 


1 


A. valida KB 


6 


6 


15 


0 


4 


4 


5 


2 


A. austera MA 


4 


4 


4 


0 


1 


0 


1 


0 


A. austera AA 


5 


4 


4 


0 


2 


0 


1 


0 


A. austera ML 


7 


6 


3 


0 


3 


0 


6 


6 


A. millepora PI 


3 


4 


4 


16 


3 


3 


13 


2 


A. millepora NG 


5 


5 


6 


12 


7 


7 


13 


7 


A. horrida PI 


3 


8 


10 


8 


4 


2 


4 


5 


A. papillare* NG 


3 


4 


6 


9 


7 


0 


7 


5 


A. papillare* PI 


7 


3 


8 


14 


5 


3 


7 


4 


A. papillare* JP 


6 


4 


4 


4 


4 


4 


5 


1 


A. pichoni* KB 


6 


3 


8 


0 


4 


2 


3 


6 


A. pichoni* L 


6 


5 


7 


0 


2 


3 


5 


6 


A. jacquelineae* KB 


3 


3 


4 


0 


7 


3 


4 


2 


A. kimbeensis* KB 


6 


4 


10 


0 


4 


2 


7 


4 


A. tortuosa* RA 


4 


3 


4 


4 


2 


0 


2 


5 


A. kirstyae* PI 


2 


4 


7 


14 


7 


5 


5 


5 


A. spathulata* PI 


2 


3 


4 


4 


2 


0 


8 


5 


A. walindii* KB 


2 


2 


5 


3 


1 


1 


3 


1 


A. rongelapensis*RA 


6 


5 


8 


0 


3 


4 


5 


5 
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